FUEL CELL PRESSURE LOADING DURING 
HYDRAULIC RAM 



John W. Patterson 



Library 

Naval Postgrackiate School 
Monterey, Calltorn a ^ ^ 40 



fSonterey, California 





FUEL CELL PRESSURE LOADING 
DURING HYDRAULIC RAM 

by 

John W. Patterson 
June 1975 

Thesis Advisor: H. L. Power, Jr. 

Approved for public release; distribution unlimited. 

T167523 



SCCUPITY classification of this page rWhTi Dmf Enfrmd) 



REPORT DOCUMENTATION PAGE 


READ IN.STRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


«. TITLE fand Su6(ltle; 

Fuel Cell Pressure Loading During 
Hydraulic Ram 


S. TYPE OF REPORT & PERIOD COVERED 

Master's Thesis 
June 1975 


6 . PERFORMING ORG. REPORT NUMBER 


7. authorc*; 

John W. Patterson 


• . CONTRACT OR grant NUM8ERr«; 


9. PERFORMING ORGANIZATION NAME ANO ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 


to. program ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 


n. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 


12. REPORT DATE 

June 1975 


• 3. NUMBER OF PAGES 
61 


U. MONITORING AGENCY NAME & ADDRESSfl/ dl/f«rww /ram Controlling O/lleoJ 

Naval Postgraduate School 
Monterey, California 93940 


IS. security class, (cl thl» rtport) 

Unclassified 


15«. DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 



16. DISTRIBUTION ST ATEMEN T fo/ fhi# R«pofO 



Approved for public release, distribution unlimited. 



17. distribution statement (of Ihm mb^irmel mnfr^d in Block 30, U difforont from Roport) 



16. SUPPLEMENTARY NOTES 



IS. KEY WORDS (Continum on r«r«re* «id« if n»c»*»mrjr «nd Identify by block numbor) 



20. ABSTRACT (Continue on r«v«r«« «id« If noco*9mry mnd Identity by block mmtbor) 

Hydraulic ram is the phenomenon associated with the damage 
caused by high speed projectile impacts with fuel tanks. It is 
divided into five phases, entry, shock, drag, cavity and exit. 
The entry, shock and drag phases were of concern in this study. 

An analytical solution for the pressures developed in the 
shock phase was formulated based on bullet impact parameters. 
Pressures developed in the drag phase were experimentally 

DD 1 jAN^i 1473 CDiTlON OF ! NOV 6S IS OBSOLETE 
(Page 1) S/N 0102-014-6601 j 



SeCORlTV CLASSIFICATION OF THIS PACE (tFhon Dots 



UNCLASSIFIED 



ftCUKlTV CLASSIFIC«TION OF THIS 0»t. EnltfiJ 



correlated with theoretical predictions. These pressures were 
proven to be dependent on front wall wave reflections. 



DD Form 1473 
, 1 Jan 73 

S/N 0102-014-G601 



UNCLASSIFIED 



security cl ASSIFICATIOM of this PAGE(lFh#n Dsf. Enffd) 




Fuel Cell Pressure Loading 
During Hydraulic Ram 



by 



John W. ,Patterson 

Lieutenant Commander, United States Navy 
B.S., Naval Postgraduate School, 1964 
M.S.I.R., George Washington University, 1971 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 

from the 

NAVAL POSTGRADUATE SCHOOL 
June 1975 



t ibrarv 

Naval Post.^^r' - . c. 
t^onterey, u.l- ., 



ABSTRACT 

Hydraulic ram is the phenomenon associated with the dam- 
age caused by high speed projectile impacts with fuel tanks. 
It is divided into five phases, entry, shock, drag, cavity 
and exit. The entry/ shock and drag phases were of concern 
in this study. 

An analytical solution for the pressures developed in 
the shock phase was formulated based on bullet impact para- 
meters . 

Pressures developed in the drag phase were experimentally 
correlated with theoretical predictions. These pressures 
were proven to be dependent on front wall wave reflections. 
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I. INTRODUCTION 



The study of aircraft survivability has become increas- 
ingly important in recent years. This increase in impor- 
tance is due to several factors. Among these are the esca- 
lating cost and concomitant complexity of modern aircraft 
making the loss of an aircraft less acceptable. The in- 
creased complexity of modern aircraft also tends to make 
them more vulnerable to anti-aircraft weapons. At the same 
time anti-aircraft defensive capability is increasing; wit- 
ness the widespread introduction of hand held infrared 
guided missiles of the Redeye class. These factors make 
survivability analysis and design of economic and military 
importance. One of the most promising areas for survivabili- 
ty analysis is the fuel system; in particular fuel cells. 

This is because the fuel cells occupy a large percentage of 
the total volume of the aircraft and are therefore most 
susceptible to damage. Damage to the fuel cell can result 
in fuel starvation, structural damage, fires or explosion; 
any of which can result in the loss of an aircraft. 

The interaction of the projectile and the fluid in the 
fuel cell produces damage that is many times greater than 
that caused by impact with an empty cell. One of the 
mechanisms that can cause damage is fire and/or explosion 
of the fluid. Hydraulic ram refers to the complicated 
dynamic interaction of the projectile, cell walls and the 
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fuel. These interactions produce internal pressure pulses 
which load the fuel cell structural members and produce 
damage . 

Hydraulic ram is usually divided into several phases 
that are associated with different major pressure pulse 
characteristics that occur during the passage of the projec- 
tile through the fluid. The five phases of hydraulic ram 
are the entry, shock, drag, exit and cavity phases. The 
entry and exit phases occur v/hen the projectile penetrates 
the respective fuel cell walls. The shock phase is associ- 
ated v/ith the initial transfer of the projectile kinetic 
energy to the fluid in the cell. This energy transfer 
creates a strong shock wave. The drag phase occurs as the 
projectile continues through the fluid and transfers more 
energy to the fluid through viscous dissipation and pressure 
drag. The cavity phase occurs as the cavity, which was 
created in the drag phase, overexpands, is compressed, re- 
bounds and continues to oscillate until the energy is 
dissipated. The entry, shock and drag phases will always 
occur in order. The exit phase may not occur, and it may 
precede or overlap the cavity phase depending on the size of 
the tank, projectile impact velocity and other factors. 
Hydraulic ram research has been conducted by various investi- 
gators. This study was designed to extend Yurkovich's 
[Ref. 1] analytical solution for the shock phase fuel cell 
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pressure loadings and to verify experimentally the analyti 
cal solution of drag phase pressures of Lundstrom [Ref. 2] 
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II, DESCRIPTION OF HYDRAULIC RAM PHENOMENON 



In analyzing hydraulic ram all five phases must be con- 
sidered. The damage associated with hydraulic ram is depen- 
dent upon the anti-aircraft threat and the physical charac- 
teristics of the fuel cell. The threat determines the path, 
energy level and size of impacting projectiles. The fuel 
cell size, geometry, location and construction, in conjunc- 
tion with the threat, determines the dominant phase of the 
hydraulic ram phenomenon. For example, there may not be an 
exit phase for a low energy threat impacting a tank with 
sufficiently large dimensions along the projectile path. 

In any phase of hydraulic ram the damage is related to 
the energy dissipated in that phase. This makes the deter- 
mination of the energy loss during the entry phase (wall 
penetration) essential. The energy dissipated during the 
entry phase is a function of wall parameters (thickness, 
yield strength) and projectile initial energy, obliquity of 
impact, material properties, and tumbling. Studies by Soper 
[Ref. 3] have shown that the percentage of energy loss for 
untumbled, normal entries can be approximated by: 



_^i"^B _ 200 ay AT. 1 . 

lost E. m ^V.2^ 



(II-l) 



1 1 

The studies also have shown that the penetration energy loss 
is independent of projectile shape. If the projectile has 
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not tumbled and strikes normal to the wall, the damage area 
(A) approximately equals the cross sectional area of the 
projectile. 

If the projectile energy is sufficiently high, it will 
penetrate into the tank fluid and produce a strong hemis- 
pherical shock wave centered about the penetration point. 

This is the start of the shock phase. The shock phase is 
characterized by high pressures of short duration (microsecs) 
that are independent of fuel cell geometry. After the bullet 
enters the fluid, it has little effect on the shock wave 
strength or shape. Investigations have shown that shock 
phase pressures generated are a function both of the energy 
transferred to the fluid at impact and of the fluid proper- 
ties. The shock wave pressures are attenuated rapidly as 
the wave is propagated through the fluid. Therefore, the 
damage associated with the shock phase is generally limited 
to the entry wall. This damage is caused by the impulsive- 
ly generated pressure field combined with whatever damage 
(cracks, etc.) the initial penetration of the wall may have 
caused . 

An analytical formulation of the shock phase problem has 
been postulated by Yurkovich [Ref. 1] . In this formulation 
the shock formation process of the shock phase has been 
modeled as a point energy source release. The formulation 
assumes that the fuel cell v;alls are rigid, and that the 
shock radius is proportional to a constant power of time. 
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The assuiaption that the radius is proportional to a constant 
power of time has been supported by several investigators 
[Ref. 4 , 5 , 6] . 

The Rankine-Hugoniot equations are used to formulate the 
flow properties at the shock front. A fourth equation is 
provided by the Tait equation of state for isothermal com- 
pressible liquid. The fluid property changes at the shock 
front are assumed to be adiabatic. If the shock radius is 
known, a shock Mach number can be computed from energy con- 
siderations. From the shock Mach number the shock front 
pressure and density can then be computed, using the above 
equations. Reference 2 postulates the initial radius is a 
function of energy deposited given by: 

R. = ( % (II-2) 

Pc 27T 

where is the energy released in the shock phase, C is the 
speed of sound and is the density in the undisturbed 
fluid. 

The flow parameters behind the shock front are com- 
puted using the unsteady motion, one-dimensional equations, 
for conservation of mass, energy and momentum, again coupled 
with the state equation. To use these equations a power law 
density distribution was assumed. The radial velocity 
distribution was found by integrating the continuity equa- 
tion. Integration of the momentum equation then yields the 
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pressure. Integration of the energy density behind the shock 
front/ knowing the shock Mach nuniber and radius, would then 
yield the energy released during the shock phase. Since 
the energy released is assumed to be known, solving the 
energy integral for the shock Mach number is possible. Once 
the shock Mach number is known as a function of radius, the 
shock radius as a function of time can be found by integra- 
tion. It is then possible to determine the pressure, shock 
Mach number and radius, as a function of time, anywhere in 
the field. In this formulation, then, the energy released 
in the shock phase must be known. 

Since, in most cases, not all of the projectile energy 
has been released during the shock phase, another phenomenon 
of hydraulic ram becomes important. This phase is referred 
to as the drag phase. It is distinguished from the shock 
phase by the gradual acceleration of the fluid vice the 
impulsive acceleration characteristic of the shock phase. 

The energy transfer, therefore, is due to viscous dissipa- 
tion and form. drag. Drag is normally characterized by 
defining a drag coefficient. The drag coefficient of a 
projectile is not a constant but varies with velocity, 
shape and tvimbling behavior. If the projectile velocity is 
below the drag divergence Mach number, however, changes in 
velocity have only a small effect on drag coefficient size. 
The effect of shape and tumbling is then the dominant 
variable for estimating the drag coefficient for a 
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projectile. Experiments [Ref. 2] have indicated that soon 
after entry (within one foot) the projectiles v/ill begin to 
tumble. Furthermore, because of the high stagnation 
pressures on the nose of the bullet, soft nosed projectiles 
may begin to deform. This phenomenon effectively changes 
the shape of the projectile and therefore the drag coeffi- 
cient. 

The pressure disturbance in the drag phase is propagated 
through the fluid at the fluid sonic velocity. Since the 
energy transfer is more gradual, the peak pressures are 
lower during this phase. But the pulses are of much longer 
duration (hundreds of microsecs) . During this phase a cavi- 
ty is formed behind the projectile. The cavity is cylindri- 
cal in shape and its size and expansion rate are dependent 
on the rate of energy dissipation by the projectile. 

Lundstrom [Ref. 2] has formulated an analytical solution 
for the drag phase pressures. In this model a velocity de- 
cay factor B is defined as a function of the drag coefficient, 
mass, area, and density of the fluid: 

B = A/2m (II-3) 

Then the rate of energy dissipation with distance is given 
as : 

^ = m BV^ (II-4) 

dx 

where B is assumed to be independent of velocity. The flow 
field properties are then solved for by using a potential 
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function that satisfies the wave equation. A line of sources 
is distributed along the projectile path with the strength 
of the sources being determined by an energy balance. The 
effect of wall restraint on cavity growth is not considered, 
so that the analysis will not yield useful pressures during 
cavity collapse. 

The pressure wave reflections at the fuel cell boundaries 
are assumed to be similar to those of a free surface. A sys- 
tem of negative image sources is utilized which result in 
zero pressure perturbation at the walls. Pressure waves re- 
flected from the cavity surfaces are not considered. 

The next phase to be considered is the cavity oscilla- 
tion phase. This is a complex phase and the least under- 
stood. Part of the complexity is due to the many reflected 
pressure waves that have had time to be of significance. 

The cavity is formed during the drag phase as the momentum 
of the fluid becomes large enough to cause the fluid to 
separate from the projectile surface. A void is created 
which is filled with fluid vapors and air that enters through 
the wall penetration point. The cavity then grows until it 
is arrested by fluid pressures developed by the walls, and 
collapses. During the collapse the momentum of the fluid 
causes the entrapped gases and air to be compressed and the 
cavity begins to expand again. It is this alternating ex- 
pansion and contraction that produces low amplitude pressure 
pulses of long duration. The phase continues until all the 
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energy is absorbed by the fluid. High speed motion pictures 
[Ref. 7] have shown that the pressure oscillations can lead 
to eventual catastrophic failure of walls already weakened 
by the drag and exit phases. 

The exit phase will occur if the energy t-ransferred from 
the projectile during the entry, shock and drag phases has 
not dissipated all the kinetic energy of the projectile by 
the time the exit wall is reached. Pressure loadings due to 
projectile proximity to the wall, and fluid motion produced 
in the earlier phases, can cause catastrophic damage to the 
exit wall. If this occurs, the latter phase (cavity oscilla- 
tion) causes little further damage but does pump fuel from 
the tank, causing fire hazard. 

From the above discussion it is apparent that each of 
the different phases of hydraulic ram produces damage to the 
fuel tank in a distinct manner. The entry phase results in 
a hole in the tank entry V7all. The damage from the shock 
phase is usually limited to the entry wall and can lead to 
catastrophic failure. The drag phase results in lower 
pressure amplitudes, but of longer duration, which can cause 
high levels of stress on all fuel cell walls. The cavity 
phase has even lower pressure amplitude but longer pulse 
widths and is oscillating in nature. This can result in 
damage to the already weakened fuel cell walls. The exit 
phase causes damage to the exit wall, usually before the 
cavity phase. It should also be apparent that impact energy 
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and the physical characteristics of the fuel cell determine 
which phase is dominant. 
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III. HYDRAULIC RAM PRESSURE EXPERIMENT 



The Naval Postgraduate School Ballistic Range was used 
to conduct the experimental part of this study. The basic 
components of the range were a rifle, velocity measuring 
equipment, simulated fuel tank, and pressure measuring and 
recording equipment. Figure III-l is a schematic of the 
basic components and Figure III-2 is a down range view. 

The rifle was a .222 caliber Remington. Cartridges were 
loaded to obtain the 7,493 in. -lb. energy level v/ith a 45 
grain spitzer ogive nose bullet. The velocity of the bullet 
was measured by using the first three Avtron No. A914T33 
chronograph screens to ti'igger two Monsanto lOlB counters. 

Two counters were used to increase the accuracy of velocity 
measurements and for redundancy. 

The simulated fuel tank was constructed of welded .160- 
inch 7075-T6 aluminum plate on the sides, back and bottom 
walls. The front wall was .05-inch 7075-T6 aluminum. A one- 
inch hole was drilled in the front wall for bullet entry. 

The hole allowed calculation of precise impact energy. The 
top of the tank was left open. The side and bottom walls of 
the tank were welded to a frame of .5-inch thick, three-inch 
aluminum angle sections for added strength. The entry and 
exit wall were bolted to the frame. The inside dimension of 
the tank was 18xl8xl8-inches . Figure III-3 shows the tank. 

A Kistler 603A quartz pressure transducer was used to 
measure the hydraulic ram pressure pulse at several points 
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internal to the tank. The transducer was mounted in a 
hollow 1.0-inch stainless steel tube. The transducer was 
held in the tube by a nylon ring. The tube was then sealed 
with RTV. In this manner the transducer and cable were kept 
dry. The tube had a 45° angle at the tip and was mounted on 
the tank at 50° from the vertical to obtain pressure measure- 
ments along a radial line. Figure I I I- 4 shows the transducer 
and probe assembly. The probe was then secured to the top 
of the tank so that it could be positioned for measurements 
at 2, 4, 6, 8 and 10 inches along a radial line from the impact 
point. 

The charge signal from the pressure transducer was fed 
into a Kistler Model 504E charge amplifier v;ith a 545A14 
filter and displayed on an oscilloscope. A Polaroid camera 
mounted on the oscilloscope recorded the pressure trace. 

The oscilloscope was triggered by the fourth chronograph 
screen mounted one inch in front of the tank wall. 

Several firings were made at each distance for consis- 
tency and accuracy. Typical results for each distance are 
shown in Figures III-5, 6, 7. 
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Figure III-2. 



Ballistic Range. 
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Figure III-3. 



Fuel Cell and Probe Installation. 



25 





26 



Figure III-4. Disassembled Probe. 





Radius = 2 in 

Scale: 1000 psi/cm, 50 microsecs/cm 




Radius = 4 in 

Scale; 500 psi/cm, 50 microsecs/cm 



Figure III-5. Pressure vs. Time, E- = 7,493 in- lb Light 
Wall. ^ 
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Radius = 8 in 



Figure III-6. 



Pressure vs. Time, E_ — 7,493 in-lb Light 
Wall. Scale: 500 PSI/cm, 50 microsecs/cm. 
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Radius = 10 in 

Scale: 500 psi/cm, 50 microsecs/cm. 



Figure III-7. Pressure vs. Time, = 7,493 in-lb Light 

Wall. 
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IV. ANALYSIS 



The first segment of this study was to improve the com- 
puter code for shock phase pressures developed by Kappel 
[Ref. 4] . This code was based on the formulation of 
Yurkovich discussed in Section II. 

To account for the energy dissipated by the projectile 
in penetrating the entry wall the formulation of Soper was - 
added to the computer code. 

A method for computing the energy deposited during the 
shock phase was needed as previous formulations assumed 
this value was known. Studies [Ref. 4, 5] have indicated 
that the energy released in the shock phase is a constant 
percentage of the impact energy. This percentage is a func- 
tion of the impact Mach number with respect to the fluid 
speed of sound. Plots of experimentally observed changes 
in projectile velocity during the characteristic shock 
formation with projectile impact v;ere made by Power [Ref. 8] . 
This plot is reproduced as Figure IV-1. In addition, the 
corresponding shock phase energy dissipation percentage of 
the impact energy is shown. It can be seen that impact Mach 
number with respect to the tank fluid is critical in deter- 
mining the shock phase energy input. For Mach numbers less 
than sonic the input percentage is approximately 10 percent 
and for moderate supersonic speeds the energy input percent- 
age is greater than 80 percent. This provided an empirical 
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procedure for determining shock phase energy input as a func- 
tion of bullet impact parameters. 

The shock front pressures are quickly attenuated and be- 
come acoustic. Figure IV-2 shows the radius and time for 
the shock wave to become acoustic versus input energy for 
both water and fuel. To calculate pressures after the 
acoustic speed was reached, a simple expression based on 
acoustic wave theory was used. The shape of the wave was 
assumed to remain constant and the pressure amplitude varied 
as the inverse of the radius . 

Since the Yurkovich theory postulates the energy release 
as a point source, the model can be extended to an explosion 
in the fuel cell. This type of event is characteristic of 
an H.E. round impact. This was accomplished by integrating 
the energy equation over a spherical volume. In this case, 
however, the shock phase energy that is released by the explo- 
sion must be known. 

Figures IV-3, 4 and 5 show the predicted pressure dis- 
tribution for .50 caliber, 12.7mm and 14.5mm AP projectiles 
impacting JP-5 fuel. The plots illustrate the very high 
pressures generated and their rapid decay. Figures IV-6 
through 11 show the predicted shock phase pressure versus 
time at a constant radius for a .222 caliber projectile im- 
pacting water. Also shown in these figures are measured 
pressures from reference 9 and analytical drag phase 
pressures to be discussed. From these figures it can be 
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seen that the shock phase pressures could not be measured. 
This is due to the very rapid rise time of the shock phase 
pressures. The pressure transducers and associated measuring 
equipment were unable to follow the signal. This phenomenon 
was also reported by Williams [Ref. 10] . 

Experiments were undertaken to verify the drag phase 
pressures predicted by Lundstrom. Holm in reference 9 re- 
ported that the drag phase pressures could not be verified. 

In that experiment the front wall of the fuel tank was a 
heavy, .5-inch, mild steel plate. Lundstrom’ s formulation 
of the problem treats the walls as free surfaces . This is 
done by using a set of negative image sources so that the 
pressure at the wall is zero. 

The degree to which a wall may be considered rigid or 
free is characterized by the expression [Ref. 2]: 

a = p C 0/m (IV-1) 

P 

where 0 is the time constant. A rigid surface corresponds 
to a=0, and a free surface to a=°°. For a .5-inch thick 
steel plate a=1.5, while for a .05-inch thick aluminum plate 
=44.5. The heavy steel plate is then an approximation to 
a rigid surface while the .05-inch aluminum plate is an 
approximation to a free surface. 

A set of image sources that can approximate the rigid 
surface would be a positive set. Using this set of positive 
images on the front wall for the pressures reported in 
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reference 9 resulted in very close agreement. This is 
illustrated in Figures IV-6 through 11 for two energy levels 
The peak pressure versus radius is depicted in Figure IV-12. 
Close agreement is shown at medium and larger radii. The 
theory does underestimate the pressures at small radii hov;- 
ever . 

The experiment was repeated using a .05-inch aluminum 
front wall to approximate a free surface. Plots of the 
pressure distribution versus time at a constant radius are 
shown in Figures IV-13 and 14. These figures also show 
close agreement between the experiment and theory when a 
free surface is assumed. 

A plot of peak pressure versus radius compared with 
theoretical values is shown in Figure IV-15. Again, the 
theory underestimates the peak pressure for small radii. 

All of the theoretical pressures were computed using a 
drag coefficient of one. This coefficient was determined 
from the experimental data of reference 3 . The theoretical 
formulation allows a variation of drag coefficient with 
time. This would account for any tumbling of the projectile 
No measurements of the tumbling behavior were possible, so 
a constant was assumed. Photographs of the bullet 
trajectory would be required to determine tumbling distances 

The pulse shape discrepancy was due to the assumption 
that was a constant. If detailed measurements are made, 
the pulse shape will be more accurately predicted. 
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Reference 11 reports that a data base was established for the 
tumbling behavior of 14.5mra and 12.7mm projectiles. From 
these data good correlation of pressures was possible. 
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PERCENTAGE VELOCITY AND ENERGY LOSS 



A 45 grain .222 Caliber Projectile 
□ 55 grain .222 Caliber Projectile 

D 5.56 uiin Tungsten Carbide Sphere 
o 5.56 mm Aluminum Sphere 
V 3.175 mm Aluminum Sphere 
<j 1.59 mm Steel Sphere 




Figure IV-1. Shock Phase Energy Release. 
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SHOCK PHASE ENERGY INPUT - E-(in.-lb 




Radius (in.) 
Time (10 ^ sec.) 



Figure IV~2. Radius and Time for Acoustic Shock Speed. 
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50 cal Projectile 
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1.0 Radius (in.) 2.0 3.0 

FIGURE iv-3. SHOCK PHASE PRESSURE PROFILE 
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FIGURE iv-4. SHOCK PHASE PRESSURE PROFILE 
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1.0 Radius (in.) 2.0 3.0 

FIGURE iv-5. SHOCK PHASE PRESSURE PROFILE 
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FIGURE iv-6. PRESSURE COMPARISON 
Eq = 12,323 IN. LB., HEAVY WALL 
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FIGURE iv-7. PRESSURE COMPARISON 
Eo =12,323 IN. LB., HEAVY WALL 
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FIGURE iv-8. PRESSURE COMPARISON 
Eq = i2,323 IN.LB., HEAVY WALL 
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FIGURE iv-9. PRESSURE COMPARISON 
Eq = 7,493 IN.LB., HEAVY WALL 
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FI CURE iv-io. PRESSURE COMPARISON 
Eq =7,493 IN. LB., HEAVY WALL 
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FIGURE iv-ii PRESSURE COMPARISON 
Eq = 7,493 IN. LB., HEAVY WALL 
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FIGURE iv-12. PEAK PRESSURE vs. RADIUS FOR HEAVY WALL 
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FIGURE iv-13 PRESSURE COMPARISON 
Eo = 7,493 IN.LB. , LIGHT WALL 
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FIGURE iv-14 PRESSURE COMPARISON 
Eq =7,493 IN.LB., LIGHT WALL 
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FIGURE iv-15. PEAK PRESSURE vs. RADIUS FOR LIGHT WALL 



V. CONCLUSIONS 



The first purpose of this study was to extend the com- 
puter code solution for shock phase pressures. This was ac- 
complished by adding solutions for the entry phase and for 
the shock phase energy release. The attached code computes 
the shock phase pressures based on a knowledge of bullet im- 
pact parameters only. It will also compute the shock phase 
pressures for a point explosion in the fluid. 

Attempts to verify the shock phase pressures experiment- 
ally were unsuccessful. This was due to the inadequate re- 
sponse time of the pressure measuring equipment. Until 
pressure transducers with extremely rapid rise times are 
available, shock phase pressures cannot be experimentally 
determined . 

The second part of this study experimentally confirmed 
the analytical solution for the pressures produced in the 
drag phase, postulated by Lundstrom. In accomplishing this, 
the importance of entry wall reflections on tank internal 
pressures was confirmed. By using both a heavy and light 
wall in the experiment, and matching these with appropriate 
image sources in the analytical solution, the pressures were 
correlated. 

The use of a constant drag coefficient in the analytical 
solution does distort the pressure wave shape. This is 
because the use of a constant drag coefficient is not con- 
sistent with the tumbling behavior of the bullet. At 
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smaller radii the experimental pressures are higher than 
predicted. This is due to the interaction of the shock 
phase pressures with the drag phase. That is, the pressure 
transducers are partially responding to the impulsive 
pressure of the shock phase. 

The results of this study show that the first three 
phases of hydraulic ram are clearly understood and can be 
predicted analytically. This should be of some assistance 
to those responsible for increasing the survivability of 
aircraft subject to light or moderate damage. 
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CD'iPUTER PROGRAM 






THIS P^OSRAM CALCULATES THE FOLLOWING 
PART I 

1. SHOCK RADIUS VERSUS TIME 

2. SHOCK MACH NUMBER VERSUS TIME 
PART II 

1. PRESSURE VERSUS TIME AND RADIUS 
FROM IMPACT BEHIND THE SHOCK 
FOR A GIVEN FLUID AND IMPACT ENERGY. 

THIS PROGRAM ASSUMES: SHOCK RADIUS IS PROPORTIONAL TO 

TIME TO THE 0.9 POWER UNTIL ACOUSTIC, A POWER LAW 
DENSITY PROFILE BEHIND THE SHOCK, A STRONG SHOCK, 
ADIABATIC ACROSS THE SHOCK, THE UNDISTURBED FLUID 
PRESSURE IS ESSENTIALLY ZERO, SPECIFIC INTERNAL 
ENERGY IN UNDISTURBED FLUID IS ESSENTIALLY ZERO. 



C 

C 
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c 
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c 
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c 
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c 
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c 
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c 

c 
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c 
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c 

c 

c 

c 

c 

c 



SYMBOL DESIGNATION 

WB=BULLET WEIGHT (GRAINS) 

C0= SPEED 0= SOUND FOR FLUID (FPS) 

RHOO= DENSITY IN UNDISTURBED FLUID (LB/FT3) 
VI=BULLET INITIAL VELOCITY (FPS) 

DC= INCREMENT OF NON-DIMENSIONAL RADIUS, DELTA 
2HETA 

DELT= INCREMENT OF TIME (SECS) 

SIGMA= YIELD STRESS OF WALL (PSD 
R3UL=BULLET RAOIJS(IN) 

THIK=EMTRY WALL THICKNESS (IN) 

XN=EXPONENT IN STATE F.QUATION (7 FOR WATER, 10.6 
FOR FUEL) 

IFLAG=ZERU -OR A HEMISPHERICAL CASE. IFLAG=ONE 
FOR A SPHERICAL CASE. 

TMAX=MAX TIME FDR CALCULATION OF SHOCK PRESSURES 
(MIkCO SECONDS) 

E03=FMFRGY RELEASE FOR SPHERICAL EXPLOSIN (INLBS) 
EOl=IMPACT ENERGY (IN-LB) 

EO=ENERGY GIVEN UP IN SHOCK PHASE (IN-LB) 

X= PERCEtNT EO LOST AT ENTRY WALL 

XX=PERCENT OF IMPACT ENERGY LOST IN SHOCK PHASE 
VS= SHOCK MACH NO. 

ET= ENERGY/VOL. = 2-E 0/ ( R=i=* 3 ) 

Fl= NON-DIMENSIONAL PRESSURE AT THE SHOCK FRONT. 
SHOCK MACH \0. 

OODM= DERIVATIVE OF NON-DIMENSIONAL PARTICLE 
VELOCITY W.R.T. SHOCK MACH NO. VS SHOCK 
MACH NO. 

XXM=INITIAL MACH NUMBER OF BULLET WITH RESPECT TO 
fluid (V I/CO ) 

DIMENSION P5I 1 (2D) , Fl( 20) ,VS( 20) , 000M(20) 

i,RXMS(20) ,DRS(20) ,AR.M5( 20) ,DT(20) ,T(20) ,RE(20) , 

1RV( 20) ,RDOD( 20),XMS{20),RS(20),RPSI(20),RF(20) 
DIMENSION TT(30U) ,XM(300 ) , R( 300 ) , F2 ( 300 ) , PS I 2 ( 300 ) , 
1D0DM2(3 00) ,3(300) ,C( i 00) ,F( 2,20) ,P( 3,2 0) 

DIMENSION CO;NRS( 3 0 0 ? , SHOC V ( 300) ,T0F0R(300) 

DIMENSION F0R(20) ,XX( 100) ,XXM( lOO) 

COMMON BQ 
DO 500 J=1 ,20 
FORI J ) =0. 0 
500 CONTINUE 
N = 19 

READ (5,100) W3,C0,RH00,VI ,DC,0ELT,SI3MA,RBUL ,THIK,XN 
READ (5,203) I F L AG , TMAX , E03 
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READ (5,101) (VS( I ) , 1=1 ,N ) 

READ( 5, 101 J (XX( I ) , 1 = 1 ,32 ) 

READ( 5, 101 ) (XXM( I) , I = 1 ,32) 

WRITE (6,200) rt3,C0,RHQ0, VI ,D:,DELT,SI3MA,RBUL,THIK,XN 
WRITE (6,204) ir-LAG,TMAX 
DO 15 1=1, N 

15 WRITE(6,202) VS( I ) 

DC 16 1=1,32 

16 WRITE (6, 201) XX(I),XXM(I) 

PART I 

PERCENTAGE OF INITIAL ENERGY DEPOSITED AND INITIAL 
SHOCK RADIUS 

E0I=WB’!= 1.4286*12.*VI*VI/(64.4*l0.=f==:=4. ) 

X = l. -( SIGMA=^'THIK^3.1416*RBUL*RBUL/E0I ) 

E0 = X»^E0 I 

CALCULATION OF ENERGY GIVEN UP IN SHOCK PHASE 
IF (IFLAG.GT.O) GO TO 3 
X6M=VI/C0 

XE=PI F2(XBM,XXM, 32,XX) 

E0=XE-E0 
GO TO 4 

3 E0=E03 

4 R0=(E0=’-32.2«144. /(RHOO-CO=:=CC*2.*3.1416) ) *’^0.3333 3 
INITIAL VALUE OF ENERGY/VOL. 

ET0=2.*E0/ ( R0**3 .0) 

INTEGRATE ENERGY EQUATION TO FIND SHOCK MACH #, 

DENSITY AND PRESSURE RATIO. 

CALL ENERG ( I FL AG , ETO , CO , RHOO , XN , VMS , PS I 10 , F 1 0 ) 

VS0=VMS 

VSI=0.01 

CALL DEVFl ( V SO , V S I , X N , DODMO ) 

COMPUTE INITIAL TIME FOR DISTANCE RO 
T0=R0/(C0*VS0=!-'12.0) 

ASSIGN SUBSEQUENT VALUES OF SHOCK MACH NO. AT WHICH 
SHOCK RADIUS WILL BE COMPUTED 

DO 10 1=1, N 

IF (VSO-VS(D) 1,2,2 

1 GO TO 10 

2 M = N-H-1 
XMS( 1)=VS0 
GO TO 11 

10 CONTINUE 

COMPUTE SUBSEQUENT VALUES OF SHOCK RADIUS 

11 RS(1)=R0 

DO 20 J=2,M 

XMS( J) =VS ( I ) 

VSE=XMS ( J) 

CALL SUBEN ( I FL AG , CO , RHOO , XN, V SE , ET ) 

RSI J) =( 2.=!=E0/(ET) )«*0.3 33 3 

1 = 1 + 1 

20 CONTINUE 
K = M-1 

COMPUTE DELTA TIME CORRESPONDING TO SHOCK MACH NO. 

AND SHOCK RADIUS 

DO 30 1=1 ,K 

DRS( I)= RSI I + l )-RS( I ) 

ARMS! I)=2./(XMS( I )+XMS( I + l ) ) 

30 DTI I )=UKS( I ) =:=ARMS ( I ) / ( C0*12.0) 
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T(1)=T0 
DO 40 J=2,M 

40 T( J)= T( J-1)1-DT( J-1 ) 

IF(RHO0.LT.i>0.0) WRIT6(6,104) 

I F( RH00.GT.61 .0) WRITE(6,102) 

WRITE (6,103) EO 
WRITE (6,110) 

PRINT SHOCK MftCH NO. AND SHOCK RADIUS VS TIME 
DO 50 1=1 ,M 

50 WRITE (6,111) T ( I ) ,XMS( I ) , RS( I ) 

100 FORMAT (8E10.0) 

101 FORMAT (6E10.0) 

102 FORMAT ( IH,// ,39X, • ENERGY' , 15X ,' FLUID WATER') 

103 FORMAT (29X,1E20.7) 

104 FCRMAT( IH,// ,39X, • ENERGY* ,15X, • FLUID FUEL’) 

110 FORMAT (1H,//,39X • T I ME ' , 1 7 X , • M ACH • , 1 7 X , • RADI US ' ) 

111 FORMAT (29X,3E20.7) 

200 FORMAT ( / / lOE 1 2 . 5/ / ) 

201 FORMAT (15X,2E20.7) 

202 FORMAT! 15X, IE20.7) 

203 FORMAT (I2,2E10.0) 

204 FORMAT ( // 1 OX 1 2 , 20XE 1 2 . 5/ / ) 



PART II 



INITIALIZE VARIABLES AND ASSIGN VALUFS USED WITH 
NON-DIMENSIONAL PRESSURE DISTRIBUTION FUNCTION 

BBB=T(M)*1 .E06fl .05 
L = BBB 
TT( 1 ) =0.0 
XM(1)=XMS(1) 

R ( 1)=R0 
F2( 1 )=F10 
PSI2(1)=PSI10 
D0DM2( 1 )=DOOMO 
Q(1)=3.«(PSI10-1 .0) 

B=-0.25 

1=1. 0/DC +1.05 
C( 1)=0.0 
WRITE (6,301) 

WRITE ( 6,302) TT ( 1 ) , X M( 1 ) , R ( 1 ) 

WRITE (6,303) 

DO 60 J=2,N 

60 C ( J)=C( J-1 )+DC 

COMPUTE NON-D.IMENSiONAL AND DIMENSIONAL PRESSURE FOR 
INITIAL SHOCK RADIUS. 

DO 70 J = 1 , I 

F(l,J)=F2(l)+PSI2(I)-(C(J)*«(Q(l)+2.0)-l.0)*(F2(l)/ 

1 PSI2 ( 1) - B« (F2( 1 ) + XM( 1 )=OGO.M2 ( 1 ) ) ) / ( Q( 1 ) + 2 .0 ) 

P( 1 , J ) = F( 1 , J)^RHG0-C0-C0=^XM( 1)*XM( 1) / ( 32 .2-144. ) 

70 CONTINUE 
K=1 
JJ=1 
SUM=0.0 

CALL FORCE(FOR,C,P,R, I ,JJ,K) 

DO 501 J=1 , I 
SUM=SUM+FOR( J) 

WRITE (6,304) C ( J ) , P( I , J ) , F ( 1 , J ) , FOR ( J ) 

501 CONTINUE 

TOFOR( 1 )=SUM 

ITERATION TO SOLVE PRESSURE, SHOCK RADIUS, SHOCK MACH 
NO. VS TIME 

DO 80 K=2,L 
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TT(K)=TT(K-1 )+DELT 
IF (TT(K)-TO) 21,21,22 

21 XM(K)=XM(1) 

R(K)=R( 1) 

F2(K)=F2( 1 ) 

PSI2( K) =PS I2( 1 ) 

D0DM2(K)=00DM2( 1 ) 

Q(K)= Q( 1 ) 

DO 55 J=1,I 

55 P (2, J)=P( 1 , J) 

GO TO 80 

22 X=TT(K) 

INTERPOLATION OF INPJT DATA AND PART I DATA FOR 
VARIABLES IN P^ESSUPE FUNCTION 

XM(K)=PIF2(X,T,M,XMS) 

R(K)=PIF2(X,T,M,RS) 

VS5=XM( K) 

CALL SHOCK ( XN , VSS , PS I 3 , F 1 2 ) 

PSI2(K) =PSI3 
F2(K)=F12 

CALL OEVFl ( VSS,VSI ,XN,D0DM5) 

DOOM 2 ( K )=D0DM3 
Q(K)=3.-(PSI2C<)-1.0) 

WRITE (6,301) 

WRITE (6,302) TT ( K ) , X M ( K ) , R ( K ) 

WRITE (6,303) 

DO 90 J = 1 , I 

F(2,J)=F2(K)+PSI2(K)*(C(J)**(Q(K)+2.0)-1.0)*(F2(K)/ 
1PSI2(K)-B*(F2(K) +XM(< )^D0DM2(K) ))/(Q(K)^-2.0) 

P( 2, J)=F( 2, J)-RH0 0^-C0=i=C0*XM(K)=!=XM(K)/(32.2*144. ) 

90 CONTINUE 
JJ = 2 
SUH=0.0 

CALL FORCE(FOR,C ,P,R, I , JJ,K) 

DO 502 J=1 , I 
SUM=SUM+FOR( J) 

WRITE (6, 304) C(J),P(2,J),F(2,J) ,FOR( J) 

502 CONTINUE 

TOFOR(K )=SUM 
80 CONTINUE 

ACOUSTIC EXTENSION OF THE MACH WAVE. 

KK=TMAX+1 .05 
LL=L+1 

DO 610 K=LL,KK 
TT(K) =TT(K-1 ) +DELT 
XM( K) =1 .00 

R(K) = R(K-1 ) + C0*DELT*12 .0 
WRITE (6,301) 

WRITE ( 6,302) TT ( K ) , XV.( K) , R ( K) 

WRITE (6,303) 

DO 600 J=1 , I 

P(3, J) = F(2, J)*RHOO*CO^CO/ (32.2^=144.=:'R(<) > 

600 WRITE ( 6,306) C ( J ) , ? ( 3 , J ) , F ( 2 , J ) 

610 CONTINUE 

WR1TE(6,305) 

CALL CONVER(R,XM,TCFjR,TT,L) 

91 FORMAT (3E15.5) 

92 FORMAT (1E15.5) 

301 FORMAT ( 1 H , / / ,3 6X , * TI ME • , 1 7X , • MACH* , I 7X , ’ R AOI US ’ ) 

302 FORMAT (29X,3E20.7) 

303 FORMAT (IH, //, 39X, 'R /R S', 17X,'P',17X,'F',19X, ’FORCE') 

304 FCRMAT(29X,4E20.7) 

305 FORMAT! IH,//,36X,’TIME',17X,'VEL0',17X,' RADUI S' , 1 7X, ' T 
lOFORCE' ) 

306 FORMAT (29X,3E20.7) 

END 
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FUNCTION PIF2 ( X , XL I ST , N, FL 1ST ) 

SECOND ORDER INTERPOLATION 

DIMENSION XLIST (100}, FLIST (100) 

BLIF (P,Q,R,S,T) = ( ( Q-P)-(S-T)/(R-Q)fS) 

IF (X-XLIST(N) ) 2,1,1 

1 I = N-1 
GO TO 5 

2 IFIX-XL IST( n ) A, 4, 6 

4 1 = 1 

5 K = 1 

GO TO 30 

6 K = 2 

7 DO 8 I = 1 ,N 

IF (X-XLIST( I ) } 9,9,3 

8 CONTINUE 
I=N 

91=1-1 

30 BLIFl = BLIF(X,XLIST(I) ,XLIST(I + 1) ,FLIST( I) ,FLIST(I + 1) 

I ) 

10 IF (K-1 ) 11, 11,12 

11 PIF2 = BLIFl 
RETURN 

12 IF({I+2)-N) 13,13,16 

13 IF ( ( I-l)-l) 15,14,14 

14 IF(ABS(XLIST(I-l)-X)-ASS(XLIST(I+2)-X))16,15,15 

15 L = 1+2 
GO TO 17 

16 L = I-l 

17 BLIF2 = BLIF ( X , XL I S T ( I) , XL I ST ( L ) ,FL I ST ( I) , F LI ST { L ) ) 
PIF2 = BLIF (X,XLIST(I+1) ,XLIST(L), BLIFl ,BLIF2) 

18 RETURN 
END 



SUBROUTINE FORC E ( FOR , C , P , R , I,JJ,K) 

THIS SUBROUTINE COMPUTES THE FORCE ON THE ENTRY WALL 
DUE TO SHOCK PHASE FOP A FRONT WALL WITH A 
PRE-PUNCHEO HOLE. 

RHRS= NON-DIMENSIONAL RADIUS OF PRE-PUNCHED HOLE 
S= SLOPE OF STRAIGHT LINE BETWEEN PRESSURE DATA 
POI NTS 

RX= THE INTERCEPT OF R /R S , NON-0 I MENS I ON A L RADIUS, 
WHERE PRESSURE IS ZERO 

FOR(iN)= F0RCE(L3) FDR AN ANNULUS WHOSE OUTSIDE RADIUS 
IS C(N) AND INSIDE RADIUS IS C(N-l) 



DIMENSION F0R(20 ) ,C(100) ,P( 2,20) ,R(300) 
RHRS=.375/R(K) 

IF(RHRS-1.) 100,25,25 
100 IFUJ-l ) 1,1,2 
1 DO 10 N=2, I 

S = (P( 1,N)-P( 1,N-1 ) )=^»4. 

IF<P(1,N-D) 3,4,4 

3 IF(P( 1, N) ) 5,5,6 

5 F0R(N)=0.0 
GO TO 10 

6 RX=(S*C(N-1)-P( 1,N-1) )/S 
IF(RHRS.GT.C(N) ) GO TO 43 
IF(RX .GT.RHRS) GO TO 60 
PH=S*( RHRS-C(N-1 ) ) +P( 1 ,N-l ) 

F0R(N)=(.5*(P(1 ,N) -PH )+PH) *3.1 41 6«(C(N)*C(N)-RHRS 
1*RHRS)*R(K)*R(K) 

GO TO 10 

60 FOR(N)=P( 1 ,N)*.5*3.1416*(C(N)*C(N)-RX*RX )*R(K)*R(K) 
GO TO 10 

4 IF(RHRS-C(N-1) ) 40,40,41 

40 FOR (N) = ( .5*( P( 1 , N )-P( 1 , N- 1 ) ) +P ( 1 , N- 1 ) ) * 3 . 1 4 16* ( C ( N ) 
1*C(N)-C {N-1)*C( N-1 ) ) *R(K)*R(K) 
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GO TO 10 

41 IF( RHRS-C(N) ) 42,43,43 

42 PH=S-{ RHRS-C (N-l ) M-P( 1 , N- 1 ) 

FOR (M) = (.5*{P(1,N)-P4)<- PH )«3. 141 6*(C(N) *C(N )-P.HRS 
1«RHRS)-R( K)=^R( K ) 

GO TO 10 

43 FOP (N) =0.0 
10 CONTINUE 

RETURN 

2 DO 20 N=2, I 

S=(P( 2,N) -P( 2,N- 1 ) )*4. 

IF(P(2,N-D) 13,14,14 

13 IF( P( 2, N) ) 15,15,16 

15 FORI N) =0.0 
GO TO 20 

16 RX=( S*C (N-l ) -P( 2 ,N-1 ) ) /S 
IF(RHRS.GT.C( N) ) GO TO 143 
I F(RX.GT.RHRS) GO TO loO 
PH=S*( RHRS-C(N-l ) )+P( 2,N-l ) 

FOR(N) = { .5=«=( P(2,N )-PH)+PH)*3.1416-(ClN)*C(N)-RHRS 
1*RHRSI*R(K)=.'=R(K) 

GO TO 20 

160 FOR(N ) = P( 2,N )*.5*3.1416=«=( C( N»*C(N )-RX*RX I’J'R (K)*R(K ) 
GO TO 20 

14 I F( RHRS-C( N- 1) ) 140,140,141 

140 FOR(N) = { .5=!'{ P( 2, N)-^( ?,N-1 ) )+P( 2,N-1 ) )*3.1416=i'(C (N) 
1*C(N)-C (N-1 )*C( N-1 ) ) = R(K) =i'R(K) 

GO TO 20 

141 IF(RHRS-C( M) ) 142,143.143 

142 PH = S-'f'( RHRS-C(N-1 ) ) <-P(2,N-l ) 

FOR(N) = { .5*( P( 2, N) -PH) + ?H) =i'3.1416=i'(C(N)«C(N)-RHRS 
1*RHRS)*R(K)=!-R( K) 

GO TO 20 S 

143 F0R(N)=0.0 
20 CONTINUE 
25 RETURN 

END 



SUBROUT INE COM VER ( R , XM , TOFOR , TT , L) 

C 

C THIS SUBROUTINE CONVERTS VELOCITY AND RADIUS 

C TO KM/SEC AND CM RESPECTIVELY 

C 

DIMENSION R{ 300) ,XM( 300) , TO FORI 300 ) , CONR S ( 3 00 ) , SHOC V ( 
1300) , TTI300) 

DO 510 11=1, L 

CONRS ( II )=R( I I )=!=2.54 

SHOCVI I I )=XM( II) *4.9/3.281 

WRI TE ( 6 ,520) TT ( I I ) , 5H3CV (II) , CONRS ( 1 1 ) , TOFOR ( 1 1 ) 

520 FORMAT! 29X,‘tE20. 7) 

510 CONTINUE 
RETURN 
END 



SUBROUTINE S HOCK ( XN , V S , PS I 1 , F 1 ) 

S THIS SUBROUTINE COMPUTES THE NON-DIMENSIONAL PRESSURE 

C AMD DENSITY AT THE SHOCK FRONT. 

C 

F(X)=-X + 1/(1-B*( X**XN-1 ) ) 

X1=(XN*( VS**2. ) ) **{ l/XN)-0.05 
B=l/ (XN*{ VS**2 . ) ) 

5 DO 10 1=1, 200 
Xl=Xl-0.005 
IFIXl .LT.0.0) 30 TO 20 
Y=F(X1 ) 

IFIY.LE. 0.001) 30 TO 40 
10 CONTINUE 
20 DO 30 J=l,200 
Xl=Xl+0 .0025 
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IF(X1 .GT.O.OJ GC TO 5 
Y=F( XI ) 

IF(Y.LE.0.001 ) GO TO 40 
30 CONTINUE 
40 psn = xi 

Fl = ( (PSI1**XN)-1 . )/(XN=i=(VS=^«2. ) ) 

RETURN 

END 



SUBROUTINE DEVFl (VS, '/SI ,XN,DOOM) 

THIS SUBROUTINE COMPUTES THE DERIVATIVE OF THE 
NON-DIMENSIONAL PARTICLE VELOCITY WITH RESPECT TO THE 
SHOCK MACH NUMBER. 

VSA=VS+VSI 
VSB=VS- VSI 

CALL SHOCK(XN,VSA,PSIA,FA) 
call SH0CK(XN,VSB,PSIB,F5) 

D0DM=(FA-F6)/(2.-VSI ) 

RETURN 

END 



SUBROUTINE ENERG ( I FL AG , ETO , CO , RHOO, XN , VMS , PS 1 1 0 , F 1 0 ) 

THIS SUBROUTINE COMPUTES THE INITIAL SHOCK MACH NUMBER 
BY integration OF THE ENERGY EQUATION. THE EXTERNAL 
SUBROUTINE QGIO IS A GUASSIAN 10 POINT QUADRATURE. 

EXTERNAL FCT 
COMMON BQ 
C = 2.0 

IF( IFLAG.EQ.O) C=1.0 
VSG=4.0 

CALL SHOCK ( XN , V SG , ^ S I 1 0 , F 1 0 ) 

A=3.1416--M RHO0/3 2.2)'^(CO=:=CO/144.O)*(VS3«''^2.O)«PSI 10 
1*F10=.''F10=:^2.0 
BQ=3.*P SI lO-l. 

CALL QGIO (0.0,1 .0, FCT, EP) 

ETG=A*EP=i''C 

IF(ETO-ETG) 40,100,50 
40 VSG=VSG-.l 

CALL SHOCK ( XN , V SG , PS 1 10 , F 1 0 ) 

BQ=3.*PSIiO-l. 

CALL QGIO (0.0, I .0,FCT,EP) 

ETG=A*EP*C 

IF(ETG-ETO) 60,100,40 
60 VSG=VSG+.001 

CALL SHOCK ( XN , V SG , PS I 1 0, F 1 0 ) 

A = 3. 1416--:'( RHOO/32 .2) ¥( C0-C0/144.0)*( VSG«'!'2 . 0 ) =!=PS I 10 
l*F10*F10-2.0 
BQ=3.*PSI10-l. 

CALL QGIO (0.0, 1 .0,FCT,EP) 

ETG=A*EP=«'C 

IF(ETG-ETO) 60,100,100 
50 VSG=VSG+.l 

CALL SHOCK ( XN , V SG , PS I 1 0 , F 1 0 ) 

A = 3.1416- (RHOO/3 2.2)-(CO*CO/] 44.0)*( VSG=i'’!'2 . 0 ) *PS I 10 
1*F10-''F1 0=i=2 .0 
BQ=3.*PSI10-1. 

CALL QGIO (0.0, 1 .0,FCT,EP) 

ETG=A«EP*C 

IF(ETG-ETO) 50,100,90 
90 VSG=VSG-.001 

CALL SHOCK ( XN, VSG, PSI lO.FlO) 

A = 3.141 6 * { RHOO/32 .2 ) C0=^C0 /1 44 . 0 ) «( VSG**2 . 0) «PS I 10 

1«F10*F1 0=!=2 .0 
BQ = 3.=i‘PSI 10- 1 . 

CALL QGIO (0.0,1 .0, FCT, 5P ) 

ETG=A=i'EP«C 
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IF(ETG-ETO) 100,100,90 
100 VMS=VSG 
RETURN 
END 



SUBROUTINE SU8EN ( I FLAG , CO , RHOO , XN , VSE , ET ) 

THIS SUBROUTINE COMPUTES THE ENERGY BEHIND THE SHOCK 
FRONT FOR A GIVEN MACH NO. THE EXTERNAL INTEGRATION 
SUBROUTINE QGID IS A GAUSSIAN 10 POINT QUADRATURE. 

EXTERNAL FCT 
COMMON 6Q 
C = 2.0 

IF( IFLAG.EQ.O) 0=1.0 
CALL SHOCK ( XN , VSE , PS I E , F E ) 

A=3. 141 6*(RH00/52. 2) =!-'l C0=!'C0/I44. 0)*( VSE«*2. 0)*PSI E 
1*FE=:'FE*2 .0 
BQ=3.*PSIE-1 . 

CALL QGIO (0.0, 1.0, FCT, EE) 

ET = A'-!'EE4'C 

RETURN 

END 



FUNCTION FCT(X) 

THIS FUNCTION SUBPROGRAM GENERATES FUNCTIONS FOR THE 
INTEGRATION ROUTINES. 

COMMON BQ 
FCT=( 1. +X**2. 

RETURN 

END 
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